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Abstract 

Thermal Bogoliubov transformation is an essential ingredient of the thermo field 
dynamics - the real time formalism in quantum field and many-body theories at 
finite temperatures developed by H. Umezawa and coworkers. The approach to 
study properties of hot nuclei which is based on the extension of the well-known 
Quasiparticle-Phonon Model to finite temperatures employing the TFD formalism 
is presented. A distinctive feature of the QPM-TFD combination is a possibility 
to go beyond the standard approximations like the thermal Hartree-Fock or the 
thermal RPA ones. 

Among numerous outstanding achievements by N. N. Bogoliubov there is the 
well-known Bogoliubov transformation for bosons [T] and fermions |2j. This uni- 
tary transformation played a crucial role in constructing microscopical theories of 
superfiuidity and superconductivity and till now has been an extremely useful and 
powerful tool in many branches of theoretical physics. Quite unexpectedly, a new 
version of the Bogoliubov transformation appeared in the middle of the 1970s when 
H. Umezawa and coworkers formulated the basic ideas of thermo field dynamics 
(TFD) [21 H] - a new formalism extending the quantum field and many-body theo- 
ries to finite temperatures. 

Within TFD [21 S], the thermal average of a given operator A is calculated 
as the expectation value in a specially constructed, temperature-dependent state 
|0(T)) which is termed the thermal vacuum. This expectation value is equal to the 
usual grand canonical average of A. In this sense, the thermal vacuum describes 
the thermal equilibrium of the system. To construct the state |0(T)), a formal 
doubling of the system degrees of freedom is introducecQ. In TFD, a tilde conjugate 
operator A — acting in the independent Hilbert space — is associated with A, in 
accordance with properly formulated tilde conjugation rules [31 HI [5]. For a heated 
system governed by the Hamiltonian H the whole Hilbert space is spanned by the 
direct product of the eigenstates of H and those of the tilde Hamiltonian H, both 
corresponding to the same eigenvalues, i.e. H\n) = En\n) and H\n) = En\n). In the 
doubled Hilbert space, the thermal vacuum is defined as the zero-energy eigenstate 

^It is worth mentioning the general statement |B] that the effect of finite temperature can be 
included in a free field theory if one doubles the field degrees of freedom. 
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of the so-called thermal Hamiltonian T-L = H — H . Moreover, the thermal vacuum 
satisfies the thermal state condition [31 IH [3] 

A|0(T)) = ae^/2r2t|o(T)), (1) 

where a = 1 for bosonic A and a = —i for fermionic A. It is seen from ([T]) that, 
in TFD, there always exists a certain combination of A and A^ which annihilates 
the thermal vacuum. That mixing is promoted by a specific transformation called 
the thermal Bogoliubov transformation This transformation must be canonical 
in the sense that the algebra of the original system remains the same, keeping its 
dynamic. The temperature dependence comes from the transformation parameters. 

The important point is that in the doubled Hilbert space the time-translation 
operator is the thermal Hamiltonian Ti. This means that the excitations of the 
thermal system are obtained by the diagonalization oiT-L. The existence of the ther- 
mal vacuum annihilation operators provides us with a powerful method to analyze 
physical systems at finite temperatures and allows for straightforward extensions of 
different zero-temperature approximations. 

In the present note, we exemplify advantages of TFD while treating the behavior 
of atomic nuclei at finite temperatures. In particular, we will show a way of going 
beyond the thermal RPA and allowing one to treat a coupling of the basic nuclear 
modes, quasiparticles and phonons [7], at finite temperatures. This problem was 
already studied in Refs. El [10] • However, the new aspects have been revealed 
recently [IT] . 

To avoid unnecessary complications in the formulae, we consider a nuclear Hamil- 
tonian which is a simplified version of the Hamiltonian of the Quasiparticle-Phonon 
Model [7\. It consists of a mean field H^p, the BCS pairing interaction impair, and a 
separable multipole-multipole particle-hole interaction Hp^. Moreover, protons and 
neutrons are not distinguished. The Hamiltonian reads 



H — + //pair + -f/ph — 

^(Ej - A)a]„a^.„ --G 4imialiw«35W«j2m2 " 2 ^o'^^^L^Xf. (2) 

jm jimij2m2 Xfi 

where a]^ and the nucleon creation and annihilation operators, ajfn = 

{—iy^"^aj_m, and M^^ is the multipole single-particle operator of the electric type 
with multipolarity A. 

At first, we apply TFD to treat pairing correlations at finite temperature (see 
also [T^ IT3]). To this aim, we make the standard Bogoliubov m, ^-transformation 
from nucleon operators to quasiparticle operators a\ a 

^jm = ^i^jm - ^jok (^i + ^] = 1) • (3) 
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The same transformation with the same u, v coefficients is apphed to nucleonic tilde 
operators aj^, a^^, thus producing the tilde quasiparticle operators a]^ and 5^-^. 

Thermal effects appear after the thermal Bogoliubov transformation which mixes 
ordinary and tilde quasiparticle operators and produces the operators of so-called 
thermal quasiparticles 13 jm and their tilde counterparts. Following the Ojima's 

formulation of the double tilde conjugation rule for fermions (a = a) [5J we use here 
the complex form of the thermal Bogoliubov transformation: 

P]m = Xj^m + ^yjOijm {x^j + v] = 1)- (4) 

The reasons for this are given in [TT] . 

Then we express the thermal Hamiltonian in terms of thermal quasiparticle op- 
erators dlj) and require that the one-body part of the thermal BCS Hamiltonian 
^Bcs = -f^sp + -f^pair " -f^sp — -f^pair bccomes diagonal in terms of thermal quasiparti- 
cles. This yields the following expressions for uj, vf 



where Sj = {Ej — Xy + A^. The gap parameter A and the chemical potential A 
are the solutions of the equations 

^ = f E(2J + - y>j^,^ N = ^(2j + l){v]x'^ + u]y% (6) 

3 j 

where is the number of nucleons in a nucleus. 

With Uj,Vj from (|5]) the one-body part of the thermal BCS Hamiltonian reads 

jm 

One can see that the Hamiltonian Hbcs describes a system of noninteracting thermal 
quasiparticles and tilde-quasiparticles with energies Sj and —Ej, respectively. 

To determine the thermal vacuum corresponding to Hbcs; we need to fix ap- 
propriately the coefficients Xj, yj. In Refs. [121 US], the coefficients were found by 
minimizing the thermodynamic potential of the system of noninteracting Bogoliubov 
quasiparticles. Here we demand that the vacuum |0(T); qp) of thermal quasiparticles 
obey the thermal state condition ([1]) 

a,„|0(T);qp) = -te'^^^^/'^a]jO{T)-qp). (8) 

Combining ([8]) and ([3]) one gets 

-1/2 



l + exp(^)]" , x, = {l-y]f\ (9) 
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We see that the coefficients yj are the thermal Fermi-Dirac occupation factors which 
determine the average number of thermally excited Bogoliubov quasiparticles in 
the BCS thermal vacuum. Equations dS]), and are the well-known finite- 
temperature BCS equations [Hj. 

In the next stage we partially take into account the particle-hole residual inter- 
action -ffph- Now the thermal Hamiltonian reads 

^ = E ^MmP.m - ^]J,J - ^ E '^o'^ {K^X, - ^iM,^} (10) 

jm \^ 

and it can be divided into two parts - Htqrpa and Tiqph- The part Htqrpa that 
contains Hbcs and the terms with even numbers of creation and annihilation op- 
erators of thermal quasiparticles is approximately diagonalized within the Thermal 
Quasiparticle Random Phase Approximation, whereas the part "Hqph containing odd 
numbers of creation and annihilation operators is responsible for the coupling of 
TQRPA eigenvectors (thermal phonons). 

To diagonalize "Htqrpa, the following operator of thermal phonon is introduced: 

where the notation [ ]^ means the coupling of single-particle momenta ji, j2 to the 
angular momentum A with the projection /i. Now the thermal equilibrium state 
is treated as a vacuum |0(T);ph) for thermal phonons. In addition, the thermal 
phonon operators are assumed to obey bosonic commutation rules. This imposes 
some constraints on the phonon amplitudes ip,4>,f] etc. (see Ref. [11] for more 
details). 

To find eigenvalues of "Htqrpa, the variational principle is applied, i.e. we find 
the minimum of the expectation value of "Htqrpa with respect to one-phonon states 
(5^^J0(T); ph) or Ql—^\0[T);ph.) under afore-mentioned constraints on the phonon 
amplitudes. As a result we arrive at the following equation for thermal phonon 
energies Uxi: 



2A+ 1 (A).2 

(A) 

^0 nh 



y^jini ^jij^S^ ^32 1 y"3i32> ^3132^^31 1^32) 



i^J)'-^' K32. 



(-)^2 , ,2 



(11) 



where fj^j,^ is the reduced single-particle matrix element of the multipole operator 

Although at the present stage phonon amplitudes cannot be determined unam- 
biguously, the TQRPA Hamiltonian is diagonal in terms of thermal phonon operators 

'Htqrpa = E ^Ai(Ql^i<5AA*i " Q\^,iQ\^i)■ (12) 
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One can see that Htqrpa is invariant under the following thermal Bogoliubov trans- 
formation: 

with X^ - — Y^- = 1. To fix the coefficients Xxi, Yxi and finally determine the phonon 
amplitudes ip, ip, 0, 0, 77, ^, we again demand that the thermal phonon vacuum obey 
the thermal state conditioiJl. For A in ([1]), it is convenient to take the multipole 
operator M^^. Then the thermal state condition takes the form 

M,^|0(T);ph) = e^TQHPA/2T]i^tjo(T);ph). (14) 
Expressing Mx^ through phonon operators we find the coefficients Xxi, Yxi 

n.= [exp(^)-l]"^', Xx.= [l + Y^;\'^\ 

The coefficients Y^^ appear to be the thermal occupation factors of the Bose-Einstein 
statistics. Thus, the phonon amplitudes are dependent on both the types of ther- 
mal occupation numbers: quasiparticle ones (the Fermi-Dirac type) and phonon 
ones (the Bose-Einstein type). The expressions for all the phonon amplitudes 
ijj, ip, (f), 0, rj, ^ can be found in [llj. 

Once the structure of thermal phonons is determined, one can find the EX- 
transition strengths from the TQRPA thermal vacuum to one-phonon states. The 
transition strengths to non-tilde and tilde one-phonon states are related by 

$L = exp(-u;,,/T)<|.^^. (15) 

This relation is equivalent to the principle of detailed balancing connecting the 
probabilities for a probe to transfer energy u to a. heated system and to absorb 
energy u from a heated system. 

Now we are ready to go beyond TQRPA and consider the effects of the term 
Hqph which is a thermal analogue of the quasiparticle-phonon interaction [7]. It 
reads 



1 c . 

^.Ph = -2 EE;7=S7 {(^l. + ^Aj5A,.(jlj2) + (h.C.)-(t.C.)} (16) 



where the notation "(h.c.)" and "(t.c.)" stands for the items which are hermitian- 
and tilde-conjugate to the displayed ones; A/''^* is the normalization factor in the 
phonon amplitudes. The operator Bx^i{jij2) reads 



^Earlier, in we have used the other procedure to this aim. That procedure seems to be 
much less evident and more lengthy. 
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where the coefficients Xjl)^, yflj^ and ^//jj are the following: 



XV" fX\ (Y 



The term "Hqph couples states with a different number of thermal phonons. To 
take into account the phonon coupling, we consider a trial wave function of the 
following form: 



l^.(JM)) 



I M 



i -^1*1 

'"^ 0(T);ph). (17) 



_ _ Af 

Alii Alii 
A2i2 -^2*2 



It should be stressed that in (|T7|) we keep the thermal vacuum of TQRPA. It means 
that we do not consider the influence of phonon coupling on thermal occupation 
numbers. Note also that the function (fT71) contains not only non-tilde one-phonon 
components but the tilde ones as well. This is a new point in comparison with 
Ref. m]. 

The function f|T7|l has to be normalized. This demand imposes the following 
constraint on the amplitudes R, R, P, S, P: 



1 2 



i 

+ E {2 W] " + KI2 (J-)] " + 2 (^-)] ' } = 1- (18) 



Alii 
A2«2 



Since the trial function contains three different types of two-phonon components, 
there are three types of interaction matrix elements which couple a thermal one- 
phonon state with two-phonon ones 

K^:iJ^) = (0(T);ph|Q,M.^^h[gUgLji|0(T);ph), 

^'2^2^^^ = (o(r);ph|Q,M.^,p4gt^,,^Qt_^J;;jo(T);ph), 

W^^:^iJt) = (0(T);ph|Q^M.Hqp4Q^^^Q^^ji|0(T);ph). (19) 

The expressions for the matrix elements U^'J^{Ji), V^^%'{Ji), and W^^i^{Ji) via the 
phonon amplitudes ip, i^, 4> etc. can be found in pTT] . 
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Applying the variational principle to the average value of the thermal Hamilto- 
nian "Htqrpa + T^qph with respect to |\E'j,(JM)) under the normalization constraint 
(fT8|) one gets a system of linear equations for the amplitudes R, R, P, S, P. The 
system has a nontrivial solution if the energy rj^, of the state |^i,(JM)) obeys the 
following secular equation: 



det 



0, 



where 



u^{ji)ui:!^{j^ 



r\\ i 



Alii 



^Aii, 



WAin + U}X2i2 + Vv 



and 



Alii 



2 I WAiii + WA2i2 



A2i2 



+ 2( 



ai+A2+J 



^Aiii — l^A2i2 ~ 



WAiii + WAaia + Vv 



(20) 



(21) 



■ (22) 



Physical effects which can be treated with the function \^!u{JM)) and Eq. (120|) 
relate to fragmentation of basic nuclear excitations like quasiparticles and phonons, 
their spreading widths and/or more consistent description of transition strength 
distributions over a nuclear spectrum in hot nuclei. 

The authors are thankful to Dr. V. Ponomarev for valuable discussions and 
comments. 
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